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The inner parts of black-hole accretion discs shine in X-rays which can be monitored and the observed spectra can be used 
to trace strong gravitational fields in the place of emission and along paths of light rays. This paper summarizes several 
aspects of how the spectral features are influenced by relativistic effects. We focus our attention onto variable and broad 
emission lines, origin of which can be attributed to the presence of orbiting patterns - spots and spiral waves in the disc. 
We point out that the observed spectrum can determine parameters of the central black hole provided the intrinsic local 
emissivity is constrained by theoretical models. 
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1 Introduction 

During the last four decades a picture of galactic nuclei has 
emerged in which massive black holes reside in the centre. 
In active galactic nuclei (AGN) accretion discs form around 
the core and provide most of the radiation emerging from 
that region. According to this scheme, intense X-rays origi- 
nate just a few gravitational radii from the black hole hori- 
zon. Manifestation of strong gravity has been searched in 
X-ray spectra of AGN with the main aim of finding the firm 
evidence for the central black hole (BH) and determining 
its parameters, namely, mass and angular momentum (for a 
detailed exposition of the subject, see Peterson 1997; Krolik 
1999). 

Broad spectral features have been expected in X-rays on 
the basis of the model in which the iron line is formed on the 
surface of a geometrically thin, optically thick and relatively 
cold medium after irradiation by a primary source (Fabian 
et al. 1989). Nowadays, the general relativistic (GR) iron 
line profiles provide a powerful tool to measure the mass 
of the black hole in AGN as well as Galactic black-hole 
candidates. Blandford & McKee (1982) and Stella (1990) 
have developed the reverberation technique that employs a 
response of the line profile following variations of the illu- 
minating primary source. Along the same line of thought, 
Matt & Perola (1992) proposed to employ variations of the 
time-integrated line properties - equivalent width, centroid 
energy and the line width. The method was developed fur- 
ther by many authors (see Fabian et al. 2000; Reynolds & 
Nowak 2003 for reviews and references). Apart from the 
mass A/., the main parameters of these models are the func- 
tional dependency of the intrinsic emissivity over the disc 
surface, I(r), the specific angular momentum a of the black 
hole and the inclination o of the accretion disc with respect 
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to observer's line of sight. High resolution in spectral and 
time domains is crucial to accomplish the analysis. 

This paper concentrates on gravitational effects that act 
on the spectral features by smearing them and moving them 
among energy bins. In this way gravity exerts the influence 
on the ultimate form of the observed spectrum. We start with 
a brief summary of the equations describing intensity (and 
polarization) propagation in the field of a black hole. De- 
tailed exploration of broad spectral features can be consid- 
ered as a method complementary to studying the overall X- 
ray continuum (Shaffee et al. 2006) on one side and narrow 
spectral features (Turner et al. 2004) on the other extreme. 

Relativistic effects are often discussed in terms of ge- 
ometrical optics with photons travelling through the empty 
spacetime - this formulation is suited quite well to various 
flavours of 'hot spots' on the disc surface, but plasma exists 
also above the disc where it influences the observed signal 
to a certain degree. Theoretical approaches have been devel- 
oped that can tackle more complicated situations, such as 
the case of dispersive media which may be able to address, 
more accurately, simultaneous observations in mutually re- 
mote parts of the electromagnetic spectrum. 

In the second part of the paper we discuss flares and 
spots as a model for X-ray variability: multiple spots are 
created on the surface of an accretion disc following the in- 
tense irradiation. The observed signal is then modulated by 
relativistic effects. This scheme is testable and it captures 
many properties of present observations. 

Mean spectra of orbiting spots resemble those of ax- 
isymmetric rings and we briefly discuss this degeneracy. It 
is inherent to time integrated data and can be resolved by in- 
creasing sensitivity and time resolution of the observations. 

In the present context, the spots represent any kind 
of localized non-axisymmetric features residing on the 
disc surface, sharing its orbital motion and modulating the 
signal. They are thought to result from the illumination of 
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a relatively cold gas of the disc after magnetic reconnec- 
tion events. Such a possibility was proposed for the X-ray 
emitting region by Galeev et al. (1979) and developed in 
many papers (Poutanen & Fabian 1999; Merloni & Fabian 
2001; Zycki 2002; Czerny et al. 2004, and references cited 
therein). 

The flare/spot scenario can be considered as a general- 
ized model of the phenomenological lamp-post geometry 
where a point-like X-ray source is located at a specified 
height on the disc axis (Henri & Pelletier 1991; Martocchia 
& Matt 1996). Martocchia et al. (2000) and Dovciak et al. 
(2004b) demonstrated by confidence contour analysis in the 
Kerr metric that this scheme can be used to determine free 
parameters if the resolution of our observations is sufficient 
and physical mechanisms underlying the line features and 
the continuum are modeled with sufficient precision, based 
on a sound understanding of the underlying physics. 

New steps towards a detailed physical analysis of the ra- 
diation mechanisms forming X-ray spectra have been pur- 
sued by several groups (Ballantyne et al. 2001; Nayakshin 
et al. 2000; Done & Nayakshin 2001; Rozahska et al. 2002; 
Collin et al. 2004; Kallman et al. 2004; Ross & Fabian 2005) 
and adapted specifically for the orbiting spot model (Czerny 
et al. 2004; Goosmann et al. 2006). A similar line of reason- 
ing stressing that the intrinsic emission of the disc medium 
has to be accurately modelled has been now adopted also 
by Brenneman & Reynolds (2006). Although the implicit 
assumptions of the codes and numerical approaches of dif- 
ferent authors are quite divers, useful comparisons were per- 
formed within the overlapping range of the parameter space 
and can be found in the above-mentioned papers. 

More complex geometries of the emitting region should 
be also explored, such as spiral waves propagating across 
the disc. This comes in accord with recognition of the im- 
portance magnetic torques may have (Krolik et al. 2005), 
but we can safely conclude that strong gravity of the central 
BH is very likely the main agent shaping the overall form of 
the X-ray spectral features from the inner disc. 

2 Preliminaries 

AGN variability time-scales extend down to a fraction of 
an hour and even less; this is comparable with the Keple- 
rian orbital period near a massive BH. The period of mat- 
ter revolving along r = const circular trajectory is i or b = 
310(ri + a)Mf seconds as measured by a distant observer, 
where M7 is the BH mass in units of 10 solar masses. 1 

1 For the gravitational field we assume the Kerr spacetime which is de- 
scribed by metric (1) below in the text. Lengths are expressed in units of the 
gravitational radius, r g =GM./c 2 =1.48 X 10 12 M7 cm. The dimension- 
less angular momentum a adopts values in the range — 1 < a < 1. Positive 
values correspond to co-rotating motion, while negative values describe 
counter-rotation (many papers assume that the accretion disc co-rotates, 
although such an assumption may not be necessarily true). Circular orbits 
of free particles are possible above the marginally stable orbit r = r ma (a) 
for the corresponding angular momentum (Bardeen et al. 1972). 



The orbital period t or b is not much longer than the light- 
crossing time, t c , because at distances of the order of a 
few r g the bulk speed of accreted material is comparable 
with the speed of light. Other time-scales relevant for black- 
hole accretion discs - thermal, sound-crossing, and viscous 
- are typically longer than t c . Because gravitation governs 
the motion near the horizon, characteristic time-scales of 
the observed variability can be scaled with M # . There are, 
however, various subtleties; the light travel time obviously 
depends on the inclination of the system with respect to an 
observer and also on scattering events that photons may ex- 
perience in the disc corona. 

The reprocessed radiation reaches the observer from dif- 
ferent regions of the system. Furthermore, if strong-gravity 
plays a crucial role photons may even follow multiple sep- 
arate paths, joining each other at the observer. Individual 
rays experience unequal time lags - for purely geometrical 
reasons and for relativistic time dilation. The final signal is 
smeared by GR effects that become dominant in the close 
vicinity of the hole and are relevant for the source spectrum 
and its variability (Laor 1991). 

The irradiation of the disc from a fluctuating source of 
primary X-rays does not lead to instantaneous response; 
time delays occur depending on the geometrical and phys- 
ical state of the gaseous material (Blandford and McKee 
1982; Taylor 1996). A single variability event is therefore 
made of a primary flare and a complex response due to 
reflection by the disc matter at various distances from 
the hole. Here we will neglect the possibility of non- 
gravitational delays, neither we will touch an interesting 
possibility that a part of reprocessed radiation is actually 
due to fast particles impinging on the disc medium (Ballan- 
tyne & Fabian 2003; Antonicci & Gomez de Castro 2005). 
Instead we will concentrate ourselves on geometrical effects 
connected with GR. Time-resolved spectroscopy offers the 
most accurate and practical method of revealing GR effects 
which is accessible with present-day technology. 

In the case of persisting patterns (t 3> t c ) a substantial 
contribution to the variability is caused by the orbital mo- 
tion (Abramowicz et al. 1991; Mangalam & Wiita 1993). 
Two types of such patterns have been examined in more 
detail, assuming the disc obeys a strictly planar geometry: 
spots (which could be identified with vortices in gaseous 
discs; cf. Abramowicz et al. 1992; Adams, Watkins 1995; 
Karas 1997), and spiral waves (Tagger et al. 1990; San- 
buichi et al. 1994). Though it may not be feasible with 
current technology, these basic possibilities are distinguish- 
able. Namely, the difference between short-living flares 
(which die on the dynamical time-scale) in contrast to the 
enduring features can be identified in time-resolved spectra 
(Lawrence & Papadakis 1993; Kawaguchi et al. 2000), but 
in order to resolve it the geometry of the system has to be 
constrained. This is not an easy task because obscuration 
e.g. by a warped disc or off-equatorial clouds may compli- 
cate the analysis (Fukue 1987; Abrassart & Czerny 2000; 
Karas et al. 1992, 2000; Hartnoll & Blackman 2000). 
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Fig. 1 Describing the GR effects acting on radiation from a source near a Kerr black hole (a = 0.998). Photons emerge from different 
points of the equatorial plane and proceed to a distant observer (located towards top of the figure), where they are collected in a detector. 
The light experiences an energy change and focusing effects which influence the observed spectrum. Time delays between photons 
are important for time-resolved spectroscopy. Left panel: the relativistic frame-dragging is illustrated here by plotting the points of 
intersection of the rays emerging from the equatorial plane. A polar grid was used, uniformly spaced in the observer plane (inclination 
9 = 60 deg; o — 90 deg corresponds to edge-on observation, i.e. in the plain of the disc). The points of intersection form ellipses 
which exhibit progressive distortion near the hole. The frame-dragging is clearly visible (the hole rotates counter-clockwise; the radii 
of horizon and of the co-rotating innermost stable orbit are indicated by circles around the origin). Right panel: Four sets of contours 
are shown. They correspond to (i) redshift function (solid lines); (ii) time delay (dashed); (iii) emission angle with respect to the local 
normal direction (dot-dashed); and (iv) radiation flux enhancement due to gravitational lensing (different levels of shading). Knowledge 
of these quantities is sufficient to compute the GR effects that are expected in observed spectra. Both plots have been constructed in 
Boyer-Lindquist coordinates; lengths are expressed in units of GM./c 2 . 



3 GR effects on light from accretion discs 

The role of the central gravitational field is illustrated in 
figure 1. Here, a view of the disc plane is shown as seen 
from above, along the BH rotation axis (it is assumed that 
the disc co-rotates with the hole around their common axis). 
What we see is a complicated interplay of geometrical ef- 
fects originating from the high curvature and rotation of 
the BH spacetime, and the abberation effects which result 
from the orbital motion of matter. All these need to be taken 
into account when determining the expected lightcurves 
(figure 2). Proceeed to the contribution by G. Matt in this 
volume for more details on GR effects seen in X-ray line 
spectra . 

The geometrical optics approximation is adequate and 
hence the task is reduced to integration of null geodesies. 
Many people have tackled this problem by following Synge 
(1967) and assuming a fixed geometry of the spacetime. 
Such an approach is well substantiated because it is the 
black hole that is of particular interest, while the disc gravity 
plays only a secondary role on the spacetime metric (it may 
be important, however, in gamma-ray burst models in which 
a companion neutron star becomes tidally disrupted). This 
situation is useful also for purposes of illustration (the case 
of a heavy disc poses a mathematically and computationally 
more challenging task because it requires to search for a si- 
multaneous solution of Einstein's equations; see Karas et al. 
1995; Usui et al. 1998). 

The basic properties of BH X-ray spectra from the in- 
nermost regions of an accretion disc are well-known (Fabian 



et al. 2000; Reynolds & Nowak 2003). The so called disc- 
line scheme has been widely applied because it allows us to 
study how the spectral characteristics are formed when ac- 
cretion proceeds in strong gravity (Fabian et al. 1989). The 
main spectral components are the X-ray continuum - pri- 
mary and reflected, and the lines - notably, the iron Ka line 
(a doublet with the rest energy of ~ 6.4 keV) and the higher 
ionisation lines (~ 6.7-6.9 keV). The line is likely formed 
within a narrow interval of radii outside r ms with the local 
intensity I{r) decreasing outwards. However, the effects of 
the BH gravity are smeared in real observations in which 
the signal comes from different, insufficiently resolved re- 
gions. Because of the integration time, the mean spectrum 
does not distinguish an orbiting source, such as a spot, from 
the entire ring. 

The gravitational field is described by the metric of Kerr 
black hole (Misner et al. 1973): 

ds 2 = -—(dt-a sin 2 9 dej)) 2 + ^ dr 2 + £ d9 2 



sin 



dt - (r 2 + a 2 ) 



(1) 

in Boyer-Lindquist spheroidal coordinates, where functions 
A(r) and £(r, 0) are known. Self-gravitation of the accreted 
gas is not taken into account in our discussion here. Among 
the properties of the metric (1) is the fact that the horizon 
occurs at A(r) = 0. Assuming sub-maximal rotation of the 
BH, \a\ < 1, the outer radius is found at the dimension-less 
r = 1 + (1 — a 2 ) 1 / 2 and it hides the singularity from a distant 
observer. As a purely GR effect, once a ^ all particles and 
photons are dragged to co-rotation with the black hole. 
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Fig. 2 Typical profiles of photometric light-curves from an orbit- 
ing spot (Karas 1996). The flux has been normalized to the max- 
imum value. Parameters of each panel: (a) r = 3, 9 a = 80 deg 
(lensing peak dominates the lightcurve near the orbital phase 0.5); 
(b) r = 44, o = 80 deg; (c) r = 44, o = 20 deg (Doppler peak 
dominates near the phase 0.75). Each curve has been normalized 
to the maximum flux. The orbital phase captures one revolution of 
the spot near a maximally rotating black hole. 



Neglecting the disc gravity is an entirely adequate 
assumption for the inner regions of AGN accretion discs 
(Novikov et al. 1973). However, the 'standard' black hole 
disc model does not seem to reproduce the required broad- 
ening of the iron line even if it is supplemented by a local 
corona above the disc. Generalized approaches have been 
discussed and the standard picture of the disc continuum 
spectrum was reconsidered to account for self-gravity (Laor 
& Netzer 1989). Karas et al. (1995) examined the gravi- 
tational effects the disc might have on the observed line 
profiles; the increase of the equivalent width was found to 
be only moderate. 

The very existence of the innermost stable orbit is due to 
general relativity. For a co-rotating equatorial disc this orbit 
has the radius (Bardeen et al. 1972) 

r ms = 3 + z 2 - [{3- z 1 ){3 + z 1 + 2z 2 )]K (2) 

where Z\ = 1 + a + a_[a + + a_], a± = (l±a)i, and 
Z2 = (3a 2 + Zf )z . The marginally stable radius spans the 
range from r ms = 1 (for a = 1, i.e. a maximally co-rotating 
BH) to r ms = 6 (for a = 0, a static BH). Black-hole rotation 
is limited to an equilibrium value by photon recapture from 
the disc (which for the standard disc model yields an equi- 
librium value of a=0.998; Thorne 1974) and by magnetic 
torques (Krolik et al. 2005). 

For a source of light obeying purely prograde Keplerian 
motion the orbital velocity is 

-l 



(3) 



v<« - A" 1 / 2 (r 2 - 2ar 1 / 2 + a 2 ) (r 3 / 2 + a) 

Here, velocity is defined with respect to a locally non- 
rotating observer at the corresponding radius. The corre- 
sponding angular velocity is Q,(r; a) — (r 3 / 2 + a)~ 1 , which 
also determines the orbital period t or b- In order to derive 
the time and frequency measured by a distant observer, one 
needs to take into account the Lorentz factor associated 
with the orbital motion, 

(r 3 / 2 + a) AV2 



By ignoring all other sources of gravitation except the 
central black hole we assume there are no secondary bodies 
influencing gravitationally the inner disc. (This may not 
be true in case of binary BHs or in the presence of global 
magnetic fields that can warp the disc away from a unique 
plane.) The radiation field is then determined by solving 
Maxwell's equations for the electromagnetic field tensor 
and its dual in a fixed curved spacetime. In the vacuum 



we write: F^-, 



0, and *F^ V - V = 0, where the aster- 



isk denotes a dual tensor. Electric field components can 
be obtained by projection onto observer's four-velocity, 
E a = F al3 up, and an electromagnetic wave is defined as 
an approximate test-field solution of the form 



F afj = 3?e u a p e 



(5) 



Two scales are introduced at this point. The phase ip{x) is 
assumed to be a rapidly varying function, while the ampli- 
tudes u a p vary slowly. This allows us to define a wave vec- 
tor, k a = ?/>,a> which is parallelly transported along the null 
geodesies, ka-p k 13 — 0, k a k a = 0. The propagation law in 
the empty space is (Anile & Breuer 1974) 



DF, 



a/3 



20F a fj — 0, 



(6) 



r i/4 Vr 3 / 2 - 3rV2 + 2a Vr 3 + a 2 r + 2a 2 



(4) 



where 9 = — -\k a - a describes the expansion of null congru- 
ences, D = ii a V a . 

Having in mind the applications to present-day X-ray 
observations, the energy shifts, gravitational lensing and 
time delays are the principal effects which originate from 
general relativity and can be tested with data that are in our 
disposal. Polarimetric information goes beyond this limit 
and has a capability of constraining parameters of our mod- 
els at a higher level. Different approaches to polarimetry 
and specific issues of X-ray polarimetry were discussed 
by various authors, namely, Cocke & Holm (1972), Anile 
et al. (1974, 1977), Madore (1974), Bicak & Hadrava 
(1975). In strong gravity, the covariant definition of basic 
polarimetric quantities is appropriate and it was devel- 
oped in various flavours: see Breuer & Ehlers (1980) and 
Anile (1989) for the definition of the polarization tensor, 
JaPjS = \{F a! jF 1 s)- A recent discussion was given in 
Broderick & Blandford (2003). 

At first, projections J a p = J a /3 7 s u 1 u 5 = {E a Ep) 
are introduced. Four parameters Sa are then given by 
S A = \{k a u a ) 2 F A , where F A (A = 0, ... 3) are again 
constructed by projecting the polarization tensor (Anile 
& Breuer 1974). These quantities satisfy the relations 
Japu 13 — 0, J a pkP = 0, lj = u a k a and can be connected 
with the traditional definition of the Stokes parameters 
(Stokes 1852; Chandrasekhar 1960). 

The normalized Stokes parameters are si = Si/ So, 
S2 = S2/ So, and s 3 = S^/So- The degree of linear and of 
circular polarization is II; = \fs\ + s\, n c = |s 3 |, and the 
total degree of polarization II = \JV^ + H 2 . 

Upon propagation through an arbitrary (curved but 
empty) space-time, the radiation flux obeys the well-known 
relation, F. dS — F, . dS . from which one can 

7 A, cm cm A, obs obs 
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Fig. 3 Isocontours of constant light-travel time St (approximately horizontal lines) and of the redshift function 1 + z (loops) are shown 
in the equatorial plane of the Kerr BH. Two panels are shown for a non-rotating BH (a = 0, left), and for a maximally rotating BH 
(a — 1, right). For a = 0, three circles are plotted with radii of the horizon, the circular photon orbit, and the marginally stable orbit. 
The three radii appear to coincide with each other in the maximally rotating case (co-rotating orbits are assumed). In the plunging region 
(shaded) the form of the isocontours depends on the assumed motion of the material, which no longer maintains Keplerian r = const 
orbits; magnetic fields likely play a major role here. All values are given in geometrized units, observer's inclination is 8 a = 50 deg. 



relate the redshift z at the point of emission with the redshift 
at a distant observer, 



1 + Z 



(7) 



(k a u a ) ohs ' A (l + z) 2 dS' 

In the Kerr metric-thin disc case, the redshift factor z 
and the local emission angle are given by 



1 + z 



r 3/2 _ 3r l/2 + 2a 



cost) = 



an 



1/2 



(8) 



r 3/2 + a _ £ 

£ and rj are constants of motion (they are connected with 
the photon ray and exist in every axially symmetric and sta- 
tionary spacetime; see figure 3). Let us remark that general- 
izations to non-Keplerian motion of the reflecting material 
were discussed by several authors: Reynolds & Begelman 
(1997) discussed radiation from the plunging region with 
matter in free-fall motion, and Dovciak et al. (2004) imple- 
mented this possibility in the KY code. Also Fukue (2004) 
examined the non-negligible radial component of velocity. 

To explore GR effects from BH accretion discs, KY 
is currently the most versatile code available publicly and 
included in the spectral fitting XSPEC package. It allows 
its user to specify time-dependent/non-axisymmetric emis- 
sivity lit; r, cp) of the disc (e.g. spiral waves), explore 
the plunging region (e.g. falling blobs), and to vary the 
black-hole angular momentum as well as the inner edge 
of the disc as free parameters. The code has also capacity 
to study GR polarization; this option cannot be currently 
exploited in practice and we need to wait for future satel- 
lites to perform the job. However, it is interesting to note 
that a complementary task can be partially accomplished 
with the help of ground based near-infrared observations 
of the Galaxy Center flares (see also Hollywood & Melia 
1997). Future X-ray satellite data will be essential to study 
GR in the neighbourhood of accreting black holes with 
outstanding precision. 

In a typical situation of an AGN accretion disc, and es- 
pecially for small radii and intermediate to large inclination 




Fig. 4 Radius r at which the maximum energy shift, 
fl+M|r-=const> reaches its peak value as a function of the incli- 
nation O . This shows how large excursions of the blue peak of the 
line can be expected. For 6 > 60 deg the maximum occurs at the 
marginally stable orbit where the inner edge of the disc is often 
assumed (details in Pechacek et al. 2005). More complete infor- 
mation about the observed width of the line can be obtained from 
graphs of the redshift factor and the gravitational lensing, as they 
vary across the disc surface (see the paper by Matt in this volume). 



angles, the line emission comes from a small fraction of 
the orbit. This implies that for observations with a moderate 
signal-to-noise ratios, only a narrow blue horn is expected 
to appear for brief time. Figure 4 shows how large variations 
of the line energy arise at different parts of the orbit. 

Eq. (7) shows that the polarization properties of the 
disc emission are modified by the photon propagation in the 
gravitational field and they may provide additional informa- 
tion about the field structure. Since the reflecting medium 
has a disc-like geometry, a substantial amount of linear 
polarization is expected because of Compton scattering. 
The set of four Stokes parameters describes polarization 
properties of the scattered light entirely. However, in order 
to compute the observable characteristics, i.e. count rates 
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of different polarization modes, one has to combine the 
reflected component with the primary continuum. The 
polarization degree of the resulting signal depends on the 
mutual proportion of the two components. 

The idea of using polarimetry to gain additional in- 
formation about compact objects is not a new one. In this 
context it was proposed by Rees (1975) that polarized 
X-rays are of high relevance. Pozdnyakov et al. (1979) 
studied spectral profiles of X-ray iron lines that result from 
multiple Compton scattering (also Angel 1969; Bonometto 
et al. 1970; Sunyaev & Titarchuk 1985; Matt 1993). The 
possibility of detecting gravitational effects on polarized 
X-rays is very attractive as it has great potential in revealing 
black holes, but the idea still awaits precise formulation as 
well as the technology usable for its practical implemen- 
tation. Temporal variations of polarization were discussed, 
in particular the case of orbiting spots (Pineault 1977; 
Connors et al. 1980; Bao et al. 1998). Further aspects of GR 
were examined by Viironen & Poutanen (2004), Dovciak et 
al. (2004a), and Horak & Karas (2006a, b), who studied the 
role of multiple images. 

4 Spots and spirals 

In their seminal papers, Cunningham & Bardeen (1972, 
1973) studied the radiation from a star revolving in the 
equatorial plane of a maximally rotating BH. They pre- 
sented a detailed analysis of periodic variations of the 
observed frequency-integrated flux that arrives at the ob- 
server. The method turns out to be a very practical one, 
even three decades later when computer capabilities have 
increased tremendously. However, the potential of their 
idea cannot be fully exploited until the theory determines 
the intrinsic emission that emerges locally from the disc. 
Neither it is well-suited for timing analyses (current status 
of this effort is summarized by Goosmann in this volume). 

Since first attempts, which were largely limited to pho- 
tometry, many people have developed various modifications 
of the original semi-analytical approach, especially for the 
spectroscopy of accretion discs: Cunningham (1975); Ger- 
bal & Pelat (1981); Asaoka (1989); Bao & Stuchlfk (1989); 
Viergutz (1993); Rauch & Blandford (1994); Zakharov 
(1994); Bromley et al. (1997); Fanton et al. (1997); Bao et 
al. (1998); Semerak et al. (1999); Martocchia et al. (2000); 
Schnittman (2005), and others. Computational challenges 
have been largely overcome thanks to a combination of 
semi-analytical approaches and sophisticated numerical 
algorithms. The analytical part of the work is largely based 
on the remarkable property of the geodesic motion in Kerr 
metric (due to Carter 1968), which is integrable and can be 
carried out in terms of elliptical integrals (de Felice et al. 
1974; Sharp 1979; Cadez & Calvani 2005). The numerical 
approach often employs pre-computed data that are stored 
for the subsequent light-ray reconstruction (Karas et al. 
1992; Dovciak et al. 2004b; Beckwith & Done 2004). 
Whichever strategy is adopted, the ray-tracing core of the 
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code has to be connected with a radiation transfer routine 
determining the intrinsic spectrum and, indeed, the under- 
standing of the transfer problem has been tremendously 
improved in recent years. 

The idea of magnetic flares irradiating the disc surface 
(Merloni & Fabian 2001) gives a promising route towards 
a more complete physical formulation of the 'bright-spot' 
model that would be based on elementary processes and 
help reducing the excessive number of degrees of freedom. 
By employing this approach, Czerny et al. (2004) computed 
the reprocessed spectra and compared the predicted vari- 
ability with MCG-6-30-15 observational data. The actual 
size of the spot is linked with the X-ray flux which origi- 
nates in flares and depends on the vertical height where the 
flares occur above the disc plane. The induced rms variabil- 
ity is a function of the energy range of the observation and 
the model parameters including the BH angular momen- 
tum and the position of the disc inner edge. Goosmann et 
al. (2006) examined different cases with the spot size and 
height ranging from a fraction of r g to several r g . It is in- 
teresting to notice that the lamp-post geometry is a limiting 
case of the flare/spot model. These results can reproduce the 
observed features, albeit the present data do not constrain all 
parameters; for example the spot size cannot be determined. 

A complementary picture assumes extended and evolv- 
ing shapes of the reflection area on the disc. Spiral shapes 
can arise from an extended spot following its decay by 
shearing motion under various kinds of instabilities oper- 
ating in the disc. Such patterns are transient but may last 
over a substantial fraction of the orbital period, sufficiently 
long to produce observable effects. Apart from AGN discs 
they have been invoked also to explain the flares in Sagittar- 
ius A* (Tagger & Melia 2006) for which individual clumps 
rather than a continuous accretion flow seem to be responsi- 
ble. Examples and useful templates of the expected spectral 
line profiles were computed (Chakrabarti & Wiita 1993; 
Karas et al. 2001; Hartnoll & Blackman 2002; Fukumura 
& Tsuruta 2004); see figure 5 where a combination of an 
instantaneous flare and a persisting spiral was assumed. 
Here, the flare quickly perishes and the spiral dominates the 
spectrum. The spiral is described by its pitch angle arctan a 
(tightly wound spirals are described by large values of the 
pitch angle) and the contrast (3 (well-defined patterns have a 
large (3). On the inner side the spiral terminates at r — r ms . 
In each panel the spectra were computed in 10 2 x 10 2 
energy-phase grid by our Kerr metric ray-tracing code. 

5 Conclusions 

The influence of strong gravity on light provides a pos- 
sibility to study black holes. However, gravity is not the 
only agent that shapes the observed spectra; emission 
mechanisms and radiation transfer in the accreting gas 
need to be understood reliably. Remarkable progress has 
been achieved on the way to model the physical processes 
which govern the form of X-ray reflection spectra emerging 
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(a) m (c) m m 01 m 




Normalized Energy (E 6.4 keV) 

Fig. 5 An examplary set of dynamical spectra of a two-arm spiral pattern with a slowly decaying flare. Different cases correspond to 
the form of the spiral, rotation parameter of the BH, and observer's viewing angle. These translate to different amplitudes and energy 
ranges spanned by observed spectra. The model parameters can be assorted in four categories: (A) BH angular momentum: a = 
(columns denoted by letters a, C, e, g), and a = 1 (b, d, f, h). (B) Observer inclination: o = 20 deg (the rows denoted by roman 
numbers i-iii), 6 Q = 50 deg (iv-vi), and o = 80 deg (vii-ix). (C) The extent of the spiral; here, the outer edge is located at: r = 7 (i, 
iv, vii), ro = 13 (ii, V, viii), and r = 20 (iii, vi, ix). (D) The actual form of the spiral structure is defined by the pitch angle arctan a and 
the contrast from the background (3: a = 3, (3 = 8 (a, b); a = V3, (3 = 8 (C, d); a = 6, (3 = 8 (e, f); a = 6, /3 = 18 (g, h). For details 
see Karas et al. (2001). 



from different regions of accretion discs. In spite of this 
advance the current description still contains some elements 
of the phenomenological approach and the emission from 
non-axisymmetric patterns cannot be computed entirely 
'from first principles'. Nonetheless, already at its present 
formulation the model of a spotted disc helps us to con- 
strain the parameters of accreting black holes. The orbital 
radius of the peak emission and the disc inclination can be 
inferred from the variations of the observed line width and 



the centroid energy. M, can be estimated by comparing 
the measured orbital period with the value expected for the 
derived radius. The model of spirals also gives the possi- 
bility of tracing the surface structures on accretion discs. 
Once these are determined the goal of measuring black hole 
rotation will be feasible with unprecedented accuracy. 
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